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Abstract. Studying transversal structure in extragalactic jets is crucial for 
understanding their physics. The Japanese led space VLBI project VSOP has 
offered arguably the best opportunity for such studies, by reaching baseline 
lengths of up to 36,000 km and resolving structures down to an angular size of 
0.3 mas at 5 GHz. VSOP observations of the jet in 0836+710 at 1.6 and 5 GHz 
have enabled tracing the radial structure of the flow on scales from 2 mas to 200 
mas and determining the wavelengths of individual oscillatory modes responsible 
for the formation of the structure observed. We conclude that these modes are 
produced by Kelvin-Helmholtz instability in a sheared relativistic flow. Our 
results point towards the stratification of the jet and the growth of different 
modes at different jet radii. We also discuss the implications of the driving 
frequency on the physics of the active nucleus of the quasar. 



1. Introduction 

Resolving the transversal (or radial) structure in extragalactic jets is a crucial 
step in our understanding of the physics of these objects. On parsec scales, 
this has become feasible only rece ntly, using space observations with the 

VSOF0 (jLobanov &: Zensusll200ll ). These observations revealed the presence of 
a double helical structure inside the jet of 3C 273, which can be attributed to a 
combination of the helical and elliptic modes of Kelvin-Helmhol tz (KH) insta- 



bility . Numerical simulations further support this interpretation (jPerucho et al 
20061 ) . We report here on further progress of this investigation. We use VSOP 
observations of the radio jet in the quasar S5 0836+710 at 1.6 and 5 GHz to 
estimate basic physical properties of the relativistic flow, focusing specifically 
on the radial profiles of its velocity and density. 

The luminous quasar S5 0836+710 at a redshift z = 2.16 ho sts a powerful 



radio jet extending u p to kiloparsec scales ( Hummel et all 19921 ). VLBI mon- 



itoring of the source (lOtterbein et all 119981 ) has yielded estimates of the bulk 



Lorentz factor 7j = 12 and the viewing angle 6j = 3° of the flow. Presence of an 
instability developing in the jet is sugges ted by the kink st ructures observed on 
milliarcsecond scales with ground VLBI fiKrichbaum et ah 1990). 

In the V SOP image of 0836+7 10 at 5 GHz, oscillations of the jet ridge line 
are observed ( Lobanov et al. 19981 ). with a wavelength of 7.7 mas. In addition 
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Figure 1. 1.6 GHz image of the jet in 0836+710. The ridge line of the jet 
is overprinted and the 100 mas mode is identified. 



to this, a 4.6 mas periodicity is foun d in the the spectral index distribution. 
Using the vortex sheet approximation ( Hardeell200Cll ). these structures have been 
identified with the helical surface and elliptic surface modes of KH instability, 
yielding estimates of basic physical properties of the flow: the Lorentz factor 
7j ~ 11, Mach number Mj ~ 6 and jet /ambient medium density ratio ^/(pj/pa) = 
0.04. High dynamic range VLBA0 observations of 0836+710 at 1. 6 GHz show the 



prese nce of an oscillation at a wavelength as long as ~ 100 mas (jLobanov et al 



l2006l ^. Fig. □ shows a 1.6 GHz image of the jet with the ridge line of the emission 
overprinted. The image also shows the identification of the 100 mas wave. This 
wavelength cannot be readily reconciled with the jet parameters determined 
from the two shorter wavelength oscillations, indicating that the flow may have 
a complex, stratified transversal structure in which emission at lower frequencies 
originates from outer layers of the flow. 

There are two basic explanations for the fact that the longest wavelength 
does not fit within the scenario described above: 1) the jet parameters derived 
may not be exact, or 2) the approximations and assumptions used in the KH 
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lin ear theory are not stri ctly valid. In this work we use the parameters derived 
Qtterbein et al. ( 19981 ) from the spectral evolution in the jet, and relax the 



m 



approximation to vortex sheet contact disconti nuity between the jet a nd the am- 

( 19981 ) ■ exploring 



bient medium, used in the stability analysis in lLobanov et al 
the possibility of the jet being sheared. 

We show that the presence of a shear layer allows to fit all the observed 
wavelengths within a single set of parameters, assuming that they are produced 
due to KH instabilities growing in a cylindrical outflow. In this picture, the 
longest mode corresponds to a surface mode growing in the outer layers, whereas 
the shorter wavelengths are identified with body modes developing in the inner 
radii of the jet. 

In Section 2, we describe the method used to solve the linear stability prob- 
lem for cylindric relativis tic flows and provide the respective solutions for the 
set of parameters given in Lobanov et al. ( 19981 ). These solutions are compared 
in Section 3 with the wavelengths observed in the jet of 0836+710. Section 4 
summarizes the main results of this work and puts them in the broader context 
of the physics of relativistic outflows. 



2. Linear Analysis 

We describe the flow in cylindrical coordinates (r, z, </)), with r and z defining the 
radial and axial directions, respectively. We consider a sheared transition in the 
axial velocity, Vz and the rest mas s density, p, between the jet and the ambient 
medium of the following form (see lPerucho et ani2007l for a detailed description 
of the linear analysis) a(r) = aoo + (ao~'2oo)/ cosh(r"'), where a(x) is the profiled 
quantity {vz or p) and oq, doo are its values at the jet symmetry plane (at r = 0) 
and at r ^ oo, respectively. The integer m controls the steepness of the shear 
layer. In the li mit m — > oo, t he configuration turns into the vortex-sheet case, 
as described bv lHardei ( 2000l ). The jet and the ambient medium are in pressure 



equilibrium. An adiabatic perturbation of the form oc (^(r) exp[i(A:2Z — wt)] 
is introduced in the equations of the Ho'rQi where a; and kz are the frequency 
and longitudinal wavenumber of the perturbation, respectively, and the function 
5f(r) defines the radial structure of the perturbation. The units used are the jet 
radius, and the speed of light, c. After some algebra, we obtain the following 
second order diff erential equation for the pressure perturbation (equation 13 in 
Birkinshawlll984l ^: 



p// ^ I 2-^lv'^^(}lz - ^VQz) _ re,U j p, ^ 
^ \ W - V^zkz " .nil 




7o 2 - wwoz) I -Pl = 0. 



Equation [T] is solved by applying the shooting method ( Press et al.lll997] ). This 
method requires integration of the equation from the jet axis to a point out- 
side the jet. We impose boundary conditions on the amplitude of the pressure 



*We assume that the magnetic field is not dynamically important. 
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Figure 2. Real (upper curves in each plot) and imaginary (lower curves in 
each plot) parts of the wavenumber versus frequency for the helical modes 
and the parameters given for the jet in 0836+710. The upper panel shows 
the solutions for m = 8 and the lower plot shows the solutions for m = 200. 
The vertical dashed lines indicate the minimum growth lengths of the surface 
modes, whereas the vertical dotted lines indicate the minimum growth length 
of the first body modes for m = 200 and the minimum growth length of all 
body modes for to = 8. The horizontal lines indicate the wave number at 
which the minimum growth length occurs. S stands for the surface mode, Bl 
for the first body mode, and B5-BQ indicate the body mode that gives the 
minimum growth lengths for to = 8. 



perturbation and on its first derivative on the jet axis, depending on the sym- 
metry or antisymmetry of the perturbation on the jet axis. Then, for a given 
frequency (uj), a complex wavenumbe r is given. The in tegration is done using 
a variable-step Runge-Kutta method ( Press et al. 19971 ). When the integration 



reaches the given point outside the jet, the value obtained is compared with the 
boundary conditions there: The Sommerfeld condition, requiring the solutions 
to approach zero at infinity, and the requirement of no incoming waves at infin- 
ity. These conditions ensure that the solutions can be given by Bessel functions 
inside the jet and Hankel functions outside. Appropriate values of the complex 
wavenumber that fit into those conditions are searched using the Miiller method 
( Press et al.l Il997l ). As stated above, the solutions are obtained in the spatial 



domain, i.e., assuming tu is real and k is complex. In this description, the inverse 
of the imaginary part of the wavenumber, 1/I'!^{k), gives the growth length (or 
e- folding length), defined as the distance at which the perturbation amplitude 
increases by on e expo nential factor. We ad opt th e jet parameters derived in 
Otterbein et~all (|l998l ): iLobanov et al.l (|l998l . l2006l ) and solve the pressure per- 



turbation equation with several different values of m. The resulting solutions 
are plotted in Fig. [2] for the helical instability modes and for two values of m: 
m = 8 (left panels) and m = 200 (right panels). The m = 8 case corresponds 
to a broad shear layer of ~ 0.6 i?j in width, whereas the m = 200 one implies a 
shear layer width of -"^ 0.1i?j, i.e., approaching the vortex-sheet case. 

The plots in Fig. [2] indicate that, in thicker shear layers (m = 8), the 
surface mode grows faster than the body modes. The minimum growth lengths 
(maximum kim) of the surface mode, indicated by vertical dashed lines in Fig. [21 
are realized at long wavelengths (horizontal dashed lines show the wavenumber 
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of the mode): Ag jj g ~ 160 at to = 0.025 c/ Rj . The smahest growth lengths 
of all body modes are achieved by the fifth and sixth order body modes at 
Ab6,8 ~ 

For a thin shear layer (m = 200) the minimum growth lengths of low order 
helical body modes are similar to that of the helical surface mode, with the min- 
imum realized for the first body mode at Abi,200 ~ 25 Ry The minimum growth 
length of the surface modes shifts to higher frequencies and shorter wavelengths 
(Agjj2oo WO Rj) compared to the respective values obtained for m = 8. For 
m < 8, the growths of all the modes are strongly reduced as the width of the 
shear layer is increased. For m > 200 the solutions converge to the vortex-sheet 
case also at the highest frequencies/shortest wavelengths. 



3. Results: From Theory to Observations 

The solutions of the linear stability problem obtained in the previous section 
can be compared with the wavelengths observed in the jet in 0836-1-710. This 
comparison requires an estimate of the radius of the jet, in order to convert the 
frequencies and wavenumbers of the solution into physical unit s. The true jet 
radius can b e estimated from the VLB A radio maps presented in iLobanov et al 



(I19981 . 120061 1. using the following relation: i?j(mas) = 0.5{Df - b^)^/'^, where A 



is the observed width of the jet and b is the beam width transversally to the 
direction of the jet. We mea sure the diamete r of the jet and the beamwidth at 
1% of the peak emission f,Wehrle et al. 19921 ) at the base of the jet, implicitly 



assuming that the outer parts of the shear layer emit less than 1% of the radio 
power at these frequencies. This yields the jet radii of -Rj,i.6 ~ 17 mas and 
-Rj,5 ~ 0.64 mas at 1.6 GHz and 5 GHz, respectively. If the jet is assumed to 
be sheared, the radius we need to consider is that at 1.6 GHz, as it includes 
the outer, slower layers of the jet. Thus, taking the radius of the jet to be 
i?j ~ 17 mas, the 100 mas structure turns into A ~ 6Rj, and the 7.7 mas into 
A ~ 0.2 Rj. 

The intrinsic (rest frame) wavelengths derived from the solutions to Eq. [1] 
shown in Fig. [2l Aint, can then be transformed to the observer's frame and ex- 
pressed also in units of Rj . The observed wavelengths are obtained by adding the 
corrections for relativistic motion, projection and cosmological effects through 
the following relation: 

\ - ^int singj 

(l + z)(l-Wccos0j)' 

where Aobs is the observed wavelength, z is the redshift, v„ is the wave speed 
and 9j is the jet angle to the line of sight. The wave speed v^j = uj/k is obtained 
from the solutions to the linear problem (Fig. [2]) . 

We focus first in the helical surface mode, most likely responsible for the 
longest observed wavelength of 100 mas. The wave speed of the helical surface 
mode ranges from ~ 0.6 c to fw ~ 0.8 c, for the thicker and thinner shear- 
layer, respectively. We obtain, in the case of the thick shear layer, a wavelength 
of 6.6 Rj, whereas in the case of a thin shear layer we obtain 8.2 Rj. 

For the case of the body modes, visible at 5 GHz, we find a wave speed of 
~ 0.2 c for the fifth and sixth body modes (those with the fastest growths) in 
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the case of a thick shear layer, and a wave speed of ~ 0.7 c for the first body 
mode in the case of a thin shear layer. These modes would produce structure 
with an observed wavelengths of 0.29 Rj and lARj, respectively. 

The wavelengths obtained for the thick shear layer are close to those ob- 
served: 6.6 Rj from theory compared to 6Rj from the observations for the 
longest structure and 0.29 Rj from theory compared to 0.2 Rj from the obser- 
vations. Thi s agreement is also re markable in the case of the 4.6 mas structure, 
identified in Lobanov et al. ( 20061 ) as an elliptic mode and not discussed here. 



This leads us to conclude that the jet is transversally stratified (i.e., has a thick 
shear layer), although the lack of error estimates precludes us from ruling out 
firmly the scenario with a thin shear layer. 



4. Discussion 

The description developed in this work explains all three major oscillations ob- 
served in the jet of 0836-1-710, if the flow speed and density are both transversely 
stratified. The longest observed wavelength of 100 mas represents the helical sur- 
face mode, whereas the two shorter ones (7.7 mas and 4.6 mas) are identified with 
the fifth/sixth order helical and elliptic body modes, respectively. The former 
is observed in the images of the jet at 1.6 GHz (see Fig. [1]) and the latter are 
observed at 5 GHz. The whole picture fits into the theory if the radius of the jet 
at the lowest observed frequency is taken, implying that the body modes have 
their maxima at inner regions of the jet than the surface mode. If we take the 
radius of the jet at 5 GHz, the short observed wavelengths may be interpreted 
as surface modes and the longes t mode cannot, in that case, be explained by 
KH theory ( Lobanov et al. 20061 ). A second implication of this result is that 



the radiating particles emitting at 5 GHz occupy a central spine in a stratified 
jet, while the longer-wavelength structure seen at 1.6 GHz is generated in outer 
layers of the flow. 

One possible alternative to the stratified jet scenari o is to conside r devi- 
ation s of the basic jet parameters from those derived in Lobanov et al.l ( 19981 . 



20061 1. To assess this alternative, we have calculated the solutions of the stability 



equation for different sets of parameters, varying the jet Lorentz factor, rest-mass 
density ratio and specific internal energy. We conclude that the jet in 0836-1-710 
must have (on axis) a Lorentz factor close to 12 or slightly smaller, a density 
ratio of 10^^-10^^ and a sound speed Cs. i ^ 0.2 or slightly large r — not deviat- 
i ng niuch fro m the parameters derived in Otterbein et al.l ( 19981 ) : lLobanov et al.l 



( 19981 . [200^ ). which supports the conclusions derived in this paper. 



Further work in the stability analysis would focus on solving the differential 
equation for the pressure perturbation (Eq. 1) using a relativistic equation 
of state that allows for inclusion of different families of particles. This will 
enable testing the solutions for a proton-electron jet and an electron-positron 
jet, providing an useful insi ght into the jet corn position. 

It has been shown by iHardee et al. ( 19941 ) that the helical surface mode 



of KH instability can be driven by an external periodic process. In this case, 
coupling this process to the helical surface mode requires the driving frequency 
to be lower than the resonant frequency. For 0836-1-710, the driving frequency 
of the helical surface mode at its minimum growth length is 0.025 c/iij would 
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imply a driving period of T^r ~ 5 .6 • 10^ yrs. This is s imilar to the driving period 
of ~ 2 • 10^ yrs found for 3C449 ( Hardee et al. 19941 ) . Periodic variations of the 



jet ejection axis can be produced by a number of physical processe s, includin, 



a bin ary black hole and a misaligned torus as the most plausible (jAppl et al 



1996l l. However, a misaligned torus can only produced a peri od of 10^ — 10^ yrs , 



depending on the accreti on disk and black-hole prop erties (|Lu k Zhoul l2005l ) 
in the case of 0836+710 (jPerucho Lobanovl tOOJU )). On the other hand, if 



we combine the value derived f or Thv with the estima te of the black hole mass 
in 0836+710 (2 ■ 10^ - lO^M©; iTavecchio et"aDl2000l ). we obtain a companion 
mass of 10^ — 10^ Mq, depending on the separation between the two black holes 
(10^^ — 10^^ cm). The latter is thus a possible explanation for the precession 
period obtained from our analysis. 

Finally, th e kiloparsec-scale str ucture of 0836+710 shows a large, decolli- 
mated feature ( Hummel et al. 19921 ). without any emission between it and the 
1.6 GHz jet. New observations with high dynamic range that have already been 
performed will give us more insight into the properties of this jet on the kilo- 
parsec scales and could show whether the cause of the decollimation on the 
kiloparsec scales is the disruption of the jet by the 100 mas structure, identified 
here clS cl helical surface mode. 
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